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Abstract 

The experimental time scale dependence of thermodynamic relations in out- 
of-equilibrium systems with aging phenomena is investigated theoretically by 
using only aging properties of the two-time correlation functions and the gen- 
eralized fluctuation-dissipation theorem (FDT). We show that there are two 
experimental time regimes characterized by different thermal properties. In 
the first regime where the waiting time is much longer than the measure- 
ment time, the principle of minimum work holds even though a system is 
out of equilibrium. In the second regime where both the measurement time 
and the waiting time are long, the thermal properties are completely differ- 
ent from properties in equilibrium. For the single-correlation-scale systems 
such as p-spin spherical spin- glasses, contrary to a fundamental assumption 
of thermodynamics, the work done in an infinitely slow operation depends 
on the path of change of the external field even when the waiting time is 
infinite. On the other hand, for the multi-correlation-scale systems such as 
Sherrington-Kirkpatrick model, the work done in an infinitely slow opera- 
tion is independent of the path. Our results imply that in order to describe 
thermodynamic properties of systems with aging it is essential to consider 
the experimental time scales and history of a system as a state variable is 
necessary. 

05.70.Ln, 05.20.-y, 81.40.C, 75.50.Lk 
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I. INTRODUCTION 



Glassy systems like spin-glasses and structural glasses below the glass transition temper- 
atures are out of equilibrium even on the macroscopic time scale. Thus, the slow dynamics 
of glassy systems has been a subject of continuous interest in the past years 0. Experi- 
mentally, the anomalous dynamical behaviors are characterized by slow relaxation with long 
time tail and aging phenomena. 

Among them, aging phenomena are the most striking dynamical behaviors as follows: 
Two-time quantities like the correlation functions C(t,t') explicitly depend on the time 
elapsed after the quench t' (the waiting time). If the waiting time t' is of microscopic time 
scale, the phenomena are merely transient on the way of relaxation to equilibrium. However, 
dependence on the waiting time continues even when t' is so large that one-time quantities 
like the magnetization are asymptotically close to time- independent values 0. Since these 
phenomena mean that the dynamics is not stationary, aging is a sign showing that these 
systems are out of equilibrium even in macroscopic time scale, i.e. several days or weeks. 

These aging phenomena appear also in mean-field models of spin-glasses and do not 
disappear even in the infinite waiting time limit |||4[] . In addition, the phase space of these 
models decomposes into a large number of areas separated with infinitely high free energy 
barriers. Thus, glassy systems never reach true equilibrium and hence they are beyond the 
scope of thermodynamics and equilibrium statistical mechanics. Hence, in order to describe 
thermodynamic properties of glassy systems, out-of-equilibrium thermodynamics based on 
a dynamical description without the assumption of ergodicity is necessary. 

In investigation of such anomalous dynamical behavior of glassy systems, it was found 
that aging phenomena have some universal properties. Theoretical analyses have suggested 
that there are two time regimes characterized by different dynamical properties of the two- 
time correlation function C(t, t') and the associated linear response function R(t, t') In 
the first time regime where the time difference t — t' is short compared to t', the dynamics 
looks stationary and the usual fluctuation-dissipation theorem (FDT) holds. On the other 
hand, in the second time regime where the time difference is comparable to t', aging phe- 
nomena occur, i.e. C(t,t') depends on t' apart from dependence on t — t'. In addition, it is 
known that the usual FDT between the correlation C(t, t') and the response function R(t, t') 
should be modified in a well-defined way which involves the rescaling of the temperature 
0. The modification was found to be valid not only for mean- field models but also for 
other glassy systems: spin-glass models with finite-range interactions ||, real spin-glasses 
0, structural glasses PHTOfl and a model of phase separation [[LI]]. In addition, it is known 
for some glassy systems that the correlation function C (t, t') obeys the scaling law that it 
depends on t and i! only through the value of £(£)/£(£')) where is a system-dependent 
increasing function of time. 

Aging of the correlation function and the modification of FDT imply that properties 
of the work done by modulating an external field in an isothermal process are completely 
different from properties predicted by traditional thermodynamics. In addition, the exis- 
tence of the two time scales implies that thermodynamic properties must strongly depend 
on experimental time scales. Hence, the experimental time scale dependence of the thermo- 
dynamic properties of glassy systems should be investigated to construct out-of-equilibrium 
thermodynamics for glassy systems. 
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In order to describe our results precisely, we summarize thermodynamics for an isother- 
mal process. Thermodynamics tells that when one quasistatically changes the external field 
the work needed for the change is independent of a path of changing the external field. In 
addition, the quasistatic work is equal to the change of the Helmholtz free energy. When the 
process is not quasistatic, the work is larger than the quasistatic work. This fact is called 
the principle of the minimum work and is derived from the second law of thermodynamics. 

We show that the properties described above do not hold in systems with aging. More 
precisely, there are two experimental time regimes characterized by different thermal prop- 
erties. The first regime is a time-domain where the waiting time is much longer than the 
time lapse of the process. We call the time lapse the measurement time. In this regime, 
the principle of the minimum work holds even though a system is out of equilibrium. More 
precisely, when the process is not infinitely slow, the work needed for the process is larger 
than the work for an infinitely slow process. In addition, value of the work for the infinitely 
slow process depends only on the initial state and the final state and hence it can play a 
role of a free energy. 

The second regime is the experimental time regime where the length of the measurement 
time are comparable to that of the waiting time. In this regime, for the single-correlation- 
scale systems such as p-spin spherical spin-glasses the work done in an infinitely slow op- 
eration depends on the path of changing the field even when the waiting time is infinite. 
This property form a striking contrast to the consequence of traditional thermodynam- 
ics described above. On the other hand, for the multi-correlation-scale systems such as 
Sherrington-Kirkpatrick model, the work done in an infinitely slow operation is independent 
of the path. 

In Sec. 0, we describe an isothermal process considered in this paper and introduce 
the two time regimes which characterize the experimental time scales and play a significant 
role in this paper. In Sec. |Tj , we see that in the first time regime, usual thermodynamic 
relations hold even though the system is out of equilibrium. The only difference is that the 
value of the work in an infinitely slow operation W s is different from that of the change of 



the free energy calculated from equilibrium statistical mechanics. In Sec. |V|, we present 
general discussion on properties of the work W s in an infinitely slow operation in the second 
time regime and derive conditions when W s depends on the path of changing the external 
field. In Sec. [V], by using the results obtained in the previous section, we show that for 
the single-correlation-scale systems the work in an infinitely slow operation depends on the 
path of changing the external field as a consequence of aging. Possibility of observation of 
this path-dependence is also discussed. In addition, we show that for the multi-correlation- 
scale systems W a is not path-dependent. Our results are summarized in Sec. [VT], where 
implications of our results on thermodynamics of glassy systems and experimental protocols 
to observe quasiequilibrium properties are discussed. 

II. AN ISOTHERMAL PROCESS AND TWO TIME REGIMES 

We describe an isothermal process to consider thermodynamic properties when aging 
occurs. A simple way to observe aging phenomena is through the following field cooling 
process (Fig. [|); The temperature of the heat bath is decreased in a small field H to a sub- 
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FIG. 1. The isothermal process to consider thermodynamic properties when aging occurs. The 
temperature of the heat bath is decreased in a small field Hq to a sub-critical temperature at 
time 0. After a waiting time t w , the field begins to change according to a given time-dependence 
Hq + H(t). This change of the field continues for a period of At. The waiting time t w and the 
measurement time At characterize the experimental time scales. 

critical temperature at time 0. After a waiting time t w the field begins to change according 
to a given time- dependence Hq + H(t). This change of the field continues for a period of 
At. The initial and final values of H(t) are H(t = t w ) = and H(t = t w + At) = AH. We 
refer to At, the time lapse of the change of the field, as the measurement time. The waiting 
time t w and the measurement time At characterize the experimental time scales. 

When the field is so weak that the response is linear, the work W done on the sample 
during the process is given in terms of the response function R(t, t') by 

W = - f tw+At dt^M(t), (2.1) 

where the "magnetization" M(t) is given by 

M(t) = f° dt'R(t,t')H + f* dt'R(t,t')H Q + f dt'R(t,t')H(t'). (2.2) 

J — oo JO Jt w 

Since we are interested in long waiting time behavior, the contribution of the first term of 
Eq. (|2.2p is ignored. For feasibility of showing long time behavior, we rewrite Eq. (|2.2| ) by 
integration by part to express M(t) in terms of a susceptibility x(t, t') instead of the response 
function; 
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t w +At dH(t) /•* , ,dH(f) , n , nn . 
dt ^T / dt'—^ X (t,t'), (2.3) 



where the susceptibility is defined as 



X (t,t')= / dt"R(t,t")- (2.4) 
Jt' 

In order to discuss the dependence of the work on the waiting time and the measurement 
time, we rewrite the expression of the work by the transformations: t — > s = (t — t w )/At, 
H(t) — > /i(s) = H(sAt + t w )/AH. We assume that dh(s)/ds and d 2 h(s)/ds 2 are finite, in 
order to exclude the unrealistic cases where the speed and the acceleration of changing the 
field is infinity. For example, the path such as h(s) = y/s is excluded, since dh(s)/ds and 
d 2 h(s)/ds 2 are infinite at s — 0. 

Thus, the work is reduced to 



W/(AH) 2 = -H / AH f 1 ds^-x(sAt + t w , 0) 

Jo ds 



-i'ds/J^ r si ds 2 dh ^ 



o dsi Jo ds2 
x x (s 1 At + t w ,s 2 At + t w ), (2.5) 

which implies that the dependence of the work on the experimental time scales is deter- 
mined by that of the susceptibility. Since the first term of Eq. (|2.5| ) becomes a constant 
— ifox(°°; 0)/ AH in the long waiting time limit which we are interested in, we will analyze 
properties of the second term of Eq. ( |2.5|) , i.e. the work done in the zero- field cooling process 
(Ho = 0), in the rest of this chapter for simplicity; 

W/(AH) 2 = - f 1 d Sl ^- r ds 2 ^- X (s 1 At + t w , s 2 At + t w ). (2.6) 
Jo as i Jo ds 2 

In order to discuss behavior of the susceptibility when aging occurs, we recapitulate 
the long-time behavior of the response function referred to in the previous section. It is 
known that there are two time regimes characterized by different behavior of the correlation 
function C{t,f) and the FDT §; 

• At long times t and t' such that t — t'^t', the correlation function is the function of 
only the time difference t — t', i.e. the time-translational invariance (TTI) holds. In 
addition, although the sample is out of equilibrium, the usual FDT holds as 

R(t t')- 1 9C{t ~ tl (2 7) 

R{t ~ t] -k^ dt' ' (2 ' 7) 



where C(r) is defined as lim^oo C(r+t, t). Since the properties of two-time quantities 
in this time regime are the same as that in equilibrium, this time regime is called the 
quasiequilibrium regime. 
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• Whereas, at long and well-separated times such that t — t' ~ t', aging occurs, i.e. the 
two-time correlation function C(t, t') depends on t' even in the long time limit t' — > oo. 
This implies no TTI and this time regime is called the aging regime. In addition, it is 
known that in the aging regime the FDT is modified as 

~ k^T dP~> (2 ' 8) 

where the FDT violation factor X is a function which depends on t, t' only through the 
dependence of the correlation function C(t, if) Thus, a system cannot be considered 
to be in a quasiequilibrium state, since the usual FDT is strongly violated. When only 
one correlation scale exists apart from the quasiequilibrium regime, it is known that 
the FDT violation factor is a constant. 

In order to clarify the meaning of the time region of the aging regime (i — t' ~ t'), we 
give an explicit expression of aging of the correlation function as 

lim C(T + t',t')=limC[T + w{T),iig(T)}=C{ii), (2.9) 

where g(t) is a system-dependent function which characterizes the aging regime and r is 
the time difference t — t'. The correlation function C(fi) depends on the waiting time t' 
through the value of \x. Occurrence of aging means that the limiting function C{pf) takes a 
non-trivial value such that < C(fi) < q. Here, q is the dynamical E-A order parameter. 
Thus, the time region of the aging regime is the time region where if / g(t — t') is finite. Here, 
we assume that \x and g(r) are positive. 

This definition of aging is illustrated in terms of contour plot of the correlation function 
onr-i' plain (Fig. (a)). The plot when aging occurs is completely different from that when 
aging does not occur (Fig. |^ (b)). The contour lines give the system dependent function g(t) 
which characterizes the aging regime, since the contour lines are given by // = t'/g(r) in the 
contour plot. 

We show two examples of the function g(t) and the aging regime for systems like spherical 
spin-glasses and real spin-glasses. In these systems, it is known that the correlation function 
behaves as a function of £(£)/£(£')> where the scaling function is a system-dependent 
function |]|. The function g(t) is given in terms of as g _1 (£) = £ _1 [r£(i)] — t, where r is 
a constant larger than unity. The inverse function of £ (t) exists since £ (t) is a monotonically 
increasing function. For example, 

1. When £(t) = t, 

lim = lim Z±l = lim I±MI = 1±£. (2.10) 

Thus, lim t-*oc C(t + t',if) is a function of fi and hence the aging regime is the time 
regime where t'/r is finite. 
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FIG. 2. (a) The contour plot on r — t' plain of the correlation function when aging occurs; 
C(t,t') = q(t' /t) 1 (7 = 0.25), which is the correlation function in the aging regime for the 5-spin 
spherical spin-glass model. The plot when aging occurs is completely different from that when aging 
does not occur in that the contour lines are not parallel to the horizontal axis. This contour plot 
gives the system-dependent function g(t) which characterizes the aging regime as a contour line of 
the correlation function, i.e. fx = t'/g(r). The straight contour lines in this figure imply that g(r) is 
a linear function of r. (b) The contour plot when aging does not occur; C(t + 1', t') = exp[— t/tq]. 
Since the value of the correlation function does not depend on t', the contour lines are parallel to 
the horizontal axis. 



7 



2. When f (t) = exp[(t/r ) 1 - a /(l - a)](0 < a < 1), 



j«» v =exp(r; rv_a) - (2 - n) 

Thus, lim t tj°° 1/ci C(r + t', if) is a function of /x and hence the aging regime is the time 
regime when t'/ r 1 /" is finite. 

From these results, it is shown by the definition of the susceptibility Eq. Q2.4| ) that the 
behavior of the susceptibility depends on the time regimes according to the dependence of 
the behavior of the response function. In the quasiequilibrium regime, the susceptibility 
depends only on the time difference r = t — if and is given by the correlation function as 

x(r) = -^[C(r)-C(0)] (2.12) 

where C(r) = linty^oo Cir + t', if). Whereas, in the aging regime, time translational in- 
variance (TTI) does not hold and the susceptibility is given by the correlation function 

as 

1 rC(t,t) 

X (t,t>) = — dCX(C) (2.13) 

Kb± JC(t,t') 



where X is the FDT violation factor defined by Eq. (|2.8| ). It implies that the susceptibility 
is a function of the correlation function. 

Since the difference of the two arguments of the susceptibility in Eq. (|2.6|) is (s\ — s 2 ) At, 
there are also two time regimes characterized by different dependence of the work on the 
experimental time scales At and t w . 

• When the waiting time is much longer than the measurement time (At t w ), the 
susceptibility which appears in Eq. ( p.6|) obeys Eq. (|2.12| ) since (si — S2)At S2At+t„, 
holds. Hence, from Eqs. ( ^.6[ ) and ( |2.12| ) the work in this regime is written as 



In the derivation, 



2&bT ksT Jo dsi 

x ds 2 ^p^C[{ Sl - s 2 )At}. (2.14) 
Jo dso 



l d Sl d 4^r ds 2 *P- = 1/2 (2.15) 
o ds\ Jo ds 2 



for any path h(s) is used. 

On the other hand, when both the waiting time and the measurement time are long 
(At ~ t w ), the susceptibility obeys Eq. fl2.13|) since (s x — s 2 )At ~ s 2 At + t w holds 
B| . In this regime, it is seen from Eq. ( p.l3j ) that the work is a functional of the 
correlation function which shows aging. The ambiguous relation t w ~ At is described 
explicitly in Sec. [TV] with the function g(t) defined in Eq. ( |2.9| ). 
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We refer to the former the quasiequilibrium regime and the latter case as the 

aging regime without any confusion with the time regimes characterized by the behavior of 
the correlation and the response function. 



III. THE WORK IN THE QUASIEQUILIBRIUM REGIME 

In this section, we investigate the properties of the work done on the sample when the 
measurement time At is short compared to the long waiting time t w (the quasiequilibrium 
regime) . 



A. The work done in an slow process 

First, we derive the work done in an slow operation such that the measurement time At 
is so long that the correlation function relax to a time-independent value. This process is 
formulated by taking the infinite measurement time limit At — > oo. Glassy systems are out 
of equilibrium even in such an infinitely slow process. Hence, we call a process where the 
external field is changed infinitely slowly a "slow process" instead of a "quasistatic process" . 

Since for the slow process the time lapse of the process is infinity, the work done in the 
slow process in the quasiequilibrium regime is given by taking the limit At — > oo after taking 
the infinite waiting time limit of Eq. ( 2.14Q . Since the integrand, the correlation function, 



is finite, the order of limit and integration can be changed. Thus, we can write the work W s 
for the slow process as 

W S /(AH) 2 = lim lim W/(AH) 2 



2k B T ksT Jo dsi 
x [ S1 ds 2 ^f^- lim C[(si -s 2 )At}. (3.1) 

Jo d,S 2 At^oo 

Here, in order to give the expression of the work W s , we introduce the dynamical 
Edwards-Anderson (E-A) order parameter defined as 

q= lim lim C(r + t w ,t w ). (3.2) 



Using Eq. (|2.15|) and the above definition, the work W s for the slow process is given by 



^ = -^^[C(0)-g]. (3.3) 

Therefore, the work for the slow process in the quasiequilibrium regime is the difference 
of a state function, since the right hand side of Eq. (|3.3j ) is independent of the path of 
changing the field h(s) and dependent only on the thermodynamic variables (AH and T) 
and constants intrinsic to the system (C(0) and q). It is important to note that this property 
is derived by using only the FDT. 
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On the other hand, for usual systems apart from glassy systems, one may expect that 
equilibrium is usually achieved when t w is long (t w — > oo). So, thermodynamics can be 
applied and the work for the slow process is equal to the change of the Helmholtz free 
energy, AF eq , calculated by equilibrium statistical mechanics. However, we show below that 
this naive expectation fails for glassy systems. More precisely, the work for the slow process, 
Eq. (|3.3|) , is different from the change of the Helmholtz free energy calculated by equilibrium 
statistical mechanics. It is because the glassy systems are out of equilibrium even when the 
waiting time is infinite. 

Assuming that A(t) is a physical quantity coupled to the external field H(t) linearly, 
we see from statistical mechanics with the assumption of ergodicity that the isothermal 
susceptibility xt is given by 

XT = j± f ((A 2 ) eq -(A)l q ) J (3.4) 

where (• • -) eq denotes average over the Gibbs-Boltzmann distribution and (• • -)j denotes 
disorder- average. Thus, the free energy difference AF eq due to change AH of the external 
field is 

&F eq = - { -^((A 2 } eq -(A}l q }j (3.5) 

where {{A) 2 q )j is the usual E-A order parameter. Since ergodicity is broken for glassy 
systems, the phase space decomposes into many pure states. Assuming that (• ■ -) denotes 
the thermal average in local equilibrium in the pure state a and F a denotes the free energy 
at the pure state a, we see that 

(A 2 ) eq = ZaP: q (A 2 )a, (3.6) 
(A)l q = {TaP e a q {A)a) 2 . (3.7) 

where P^ q is the probability that the system is found at the pure state a in true equilibrium 
and is given by P£ q = exp(—/3F a )/Z. 

In order to compare W s , Eq. (|3.3j), with AF eq , Eq. (|3.5|), we rewrite Eq. ( |3.3p in terms 
of pure states. Since the correlation function is given by C(t,t') = ((A(t)A(t')))j and the 
local equilibrium in pure states is achieved in the long time limit, 



C(0) = lim ((A(t') 2 ))j = J2Pa m ((A\)j, (3.8) 

v— »oo 

a 

q = hm lim ((A(r + t')A{t')))j = £ Pr((A) 2 a )j (3.9) 



T^OO t'— >00 

a 



where P^ m is the probability that the system is found in pure state a at time 0. Since the 
system is out of equilibrium at time 0, P^ q ^ P™ 1 - Thus, ((A 2 ) eq ) j is not equal to C(0). In 



addition, from Eqs. (|3.7|) and (|3.9|), {{A) eq )j is not equal to q. Consequently, we conclude 
that the work W s for the slow process does not coincide with the free energy difference 
calculated by statistical mechanics with assumption of ergodicity. 
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B. The work when the measurement time is finite 

We discuss the properties of the work when the measurement time is finite. From Eqs. 
( p. 14j ) and (|3.3j), the difference between the work when the process is not slow and the work 
for the slow process is given by 

T „ TTf (AH) 2 /-i , dh( Sl ) r'i , dh(s 2 ) r ^ r/ N A . , 

W-W s = i— -f- / d Sl — ^ / ds 2 —^{C[(s t - s 2 At - q}. 3.10 
KbT Jo dsi Jo ds 2 

C[(s\ — s 2 ) At] > q when At is finite, since the quasiequilibrium regime is considered. Hence, 
the difference from the work for the slow process is positive for any path of changing the 
field, when the measurement time is finite. This implies the principle of minimum work; 

W > W„ (3.11) 

where the equality holds only when the measurement time is infinite, i.e. the slow-process 
limit. 

Our result is very similar to the consequence of thermodynamics which tells the work for 
non-quasistatic process is larger than the change of the Helmholtz free energy. However, our 
result is different from that and beyond the scope of thermodynamics since in our discussion 
the initial state and the final state of the process are out of equilibrium and the value of the 
work for the slow process is different from the value of change in the Helmholtz free energy 
derived by equilibrium statistical mechanics. 



C. Long measurement time behavior of W — W s 



The long measurement time behavior of W — W s is analyzed in this subsection. We show 
that the behavior of W — W s when the measurement time At is long but finite is determined 
by the long time behavior of the correlation function. We describe below the results for four 
types of behavior of the correlation function which include almost all types of relaxation 
of the correlation, e.g., the exponential, the power law fI5|| , the logarithmic |16[ and the 
stretched exponential relaxation JT^|. The derivations are given in Appendix [A[ 



1. When the correlation function behaves as lim 7 
given by 



,t[C{t) — q] =0, the difference is 



W-W q 



[AHf 
k B T 



ds 



dh(s) 



ds 



1 

At' 



(3.12) 



where 



dr[C(r) - q] 



(3.13) 



and W s denotes the work for the slow process. This case includes the power law 
relaxation such that C{r) ~ q + cr~ Q when a > 1 and the stretched exponential 
relaxation (exp(— ar n ),0 < n < 1) as well as the exponential relaxation 
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2. When the correlation function behaves as lim^oo r[C(r) — q] = c, the difference is 
given by 



(AH) 2 r l 

KB-/ JO 



'/- 



n (At) 

—7-^- 3 - 14 



These two results show that — W s is proportional to the inverse of the measurement 
time only when lim^oo t[C(t) — g] = 0. In these two cases, the difference from the 
work for the slow process takes the minimum value when the field increases linearly 
as h(s) = s, since we assume that dh(s)/ds and d 2 h(s)/ds 2 is finite. 

3. When the correlation function obeys the power law relaxation, such that C(t) ~ 
q + cr~ a (0 < a < 1), the difference is given by 

___ (AH) 2 /-i d/i( Sl ) 



/cbj" Jo ds\ 

x / ds 2 r 7 ^-77-^-- 3.15 

Jo ds 2 ( Sl - s 2 ) a (At) a v ; 

In this case, the difference from W s obeys the power law whose exponent is equal to 
the exponent of the correlation function. 

4. When the correlation function obeys the logarithmic relaxation C(r) ~ q + c/ln(r), 
the difference also obeys the logarithmic relaxation as 

In this case, the difference does not depend on the path h(s). 

Experimentally, these results tell how long the measurement should take and the suitable 
path of changing the external field in order to determine the value of the work for the slow 
process, i.e. the difference of a state function in quasiequilibrium regime. 



IV. THE WORK IN THE AGING REGIME; GENERAL RESULTS 

In this section, we discuss the properties of the work in the isothermal slow process when 
both the waiting time and the measurement time are long, i.e. in the aging regime, by using 
the modified FDT and aging of the correlation function introduced in Eqs. ( |2.8|) and ( p.9|) . 
Condition when the work for the slow process depends on the path changing the external 
field is obtained. By using the results obtained in this section, the path-dependence of the 
work for particular systems is discussed in the next section, Sec. 0. 
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A. Path-dependence of the work for the slow process 



As shown in Sec. in the aging regime the work is a functional of the correlation 
function which shows aging. The slow process of the aging regime is given by taking the 
long measurement time limit At — > oo with holding the relation t w ~ At, which guarantees 
that the slow-process limit is taken within the aging regime. The slow-process limit has 
always this meaning in this section. 

The order of the slow-process limit and the integration can be changed since the suscep- 
tibility x is finite. From Eq. (|2.6| ) the work W s for the slow process in the aging regime is 
given by 



o dsi Jo ds 2 
x x [ lim C( Sl At + t W} s 2 At + t w )}. (4.1) 



At— >oo 
i™~Ai 



Thus, the dependence of limAt^oo C(s\At + t w ,s 2 At + t w ) on s\ and s 2 determines the 

t w ~At 

dependence of the work for the slow process on the path of changing the field h(s); W s is 
independent of h(s) if lim At-oo C(s 1 At + t w , s 2 At + t w ) does not depend on s± and s 2 . On 

the other hand, W s is a functional of h(s) if lim At-oo C(siAt + t w , s 2 At + t w ) depends on 

t„~At 

Si and s 2 . We derive conditions when W s depends on h(s) in the rest of this section. 

In order to analyze lim At->oo C(siAt+t w , s 2 At+t w ), we clarify the meaning of t w ~ At or 

the aging regime by using g{r) in Eq. (pT9|) . By substitution t' — > s 2 At + t w , r —> (si — s 2 )At 
in Eq. (|2.9|) , we see that the following equation holds in the region except for the point 
Si = s 2 which does not contribute to the value of the work by itself; 

lim C(s 1 At + t w ,s 2 At + t w ) = CU). (4.2) 

At— »oo 
s 2 At+t w 



* g[(s 1 -s 2 )At] 

This equation implies that the aging regime is the experimental time regime where (s 2 At + 
t w )/g[(s\ — s 2 )At] is finite and the slow-process limit of the correlation function in the 
integrand in Eq. (^4.1|) is a function of lim At^oo (s 2 At + t w )/g[(s\ — s 2 )At]. 

t w ~ At 

We discuss below the aging regimes and dependence of lim At-oo C(si At + t w , s 2 At + t w ) 

t w ~At 

on Si and s 2 in four cases of the different long-time behaviors of g(t) which exhaust all 
possibilities (see Fig. |3|). 

1. When lim^oo r/g(r) = oo, t 1 j git — t') is finite only in the region if <C t — t'. In the 
region of integration of Eq. (|4.1|) except for the point s 2 = 0, 



[s 2 At + t w )/g[( Sl -s 2 )At} ^oo (4.3) 



as At — ► oo since t w > 0. An example of such g{r) is g{r) = logr. When g{r) = logr 
and s 2 7^ 0, liniA^oo^At + t w )/g[(si — s 2 )At] = oo as shown in Eq. ([4.3|). Since the 
point s 2 = does not contribute to the value of the work by itself, the experimental 
time regime where (s 2 At + t w )/g[{s\ — s 2 )At] is finite does not exist in this case. If 
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FIG. 3. The schematic illustration of the four types of the long-time behavior of the function 
g(r). For example, g{r) = logT (type 1.), g(r) = r (type 2.), g{r) = t 2 (type 3.), g(r) = expr (type 
4.) are plotted against 1/r. It is important to note that these four types exhaust all possibilities 
of long-time behavior. 



there is no other aging regime except for that characterized by such g(r) exists, only 
the quasiequilibrium regime contributes to the value of the work and hence the work 
for the slow process is independent of h(s) and is given by Eq. Q3.3|) . 

2. When g(r) = r, 

s 2 At + t w _ s 2 At + t w _ s 2 + t w J At ^ ^ 



g[(si - s 2 )At] (s% - s 2 )At si - s 2 



Thus, from Eq. (|4.2|) , lim At-oo C(s\ At+t w , s 2 At+t w ) is a function of (s 2 + fi')/(si — 

n'=tv>/At 

lim C(s\At + t w , s 2 At + t w ) 

At— *oo 

= lim C(s\At + t w , s 2 At + t w ) 

_ S2At-\-t w 
s l~ s 2" s[(s 1 -s 2 )Ai] 

= C7(^±^). (4.5) 
S1-S2 

It implies that the aging regime is the region where t w /At is finite. In this aging 
regime the slow-process limit of the correlation function in the integrand of Eq. ( |4.1|) 
depends on s± and s 2 if f(^) in Eq. ( |4.2| ) is not a constant. It means that the work 
for the slow process, W s , depends on the path of change of the external field, h(s). 

3. When lim^oo rjg{j) = and \\vcl t -+ 00 g(r) / g(sr) = g r (s) (s < 1) exists, 
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s 2 At + t w s 2 At + t w g(At) 



g[( Sl -s 2 )At] g(At) g[( Sl - s 2 )At] 

~> ^AT) 9r{si - S2) ' (4 ' 6) 

as At — > oo. An example of such g{r) is g(r) = r a {a> 1). When g{f) = r a {a > 1), 
Eq. (p|) holds; 



s 2 At + t v 



g[( Sl - s 2 )At) g(At) (s x - s 2 )« 
as At — > oo. From Eqs. ( |4.2|) and 



(4.7) 



lim C(siAt + t w , s 2 At + t v 

At— >oo 
f'=<iu/9(At) 



lim C(siAt + t w , s 2 At + t w ) 



At^OO 



/ / , soAt-ftio 

" r(si ' S2) ^Kn-s2)At] 

= C'[/iV(si-s 2 )]. (4.8) 

It implies that the aging regime is the region where t w /g(At) is finite. In this aging 
regime the slow-process limit of the correlation function in the integrand of Eq. ( |4.1| ) 
is a function of s\ and s 2 if f(fi) in Eq. (|4.2| ) is not a constant. It means that the work 
W s depends on h(s). 

4. When lim,-^^ r/g(r) = and lim^oo g(r)/g(sr) = oo (s < 1), we prove below 
that the aging regime for experimental time scales does not exist. We assume that the 
aging regime exists. It means that a function p(t) exists such that lim At^oo (s 2 At + 

/i=i«j/p(Ai) 

t w )/g[(si — s 2 )At] is finite. Then, from the condition lim^oo r/g{r) = 0, 

s 2 At + t w 
lim — . - , 

At ir g[(si-s 2 )At] 

= hm 



At^oo 



- S 2 )At] 



>p(At) 

= ^ r/ P(A l f ( 4 - 9 ) 

Ai-+oo y[( Sl - S2 )At] 

Assuming that s' x — s' 2 > si — s 2 , we see from the condition Hm r _ j . 00 ^(r)/^(sr) = 
oo (s < 1) that 

s 2 At + t w , s^At + t^ 



lim / lim 

At 1^_ " * 2 ) Af] ^ - 4) At] 

M ~p(At) M ~p(At) 

lim MAt) £[(4 - s' 2 )At] 
At^oo g[(si — s 2 )At] ^p(At) 
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nm g[(s[ - s' 2 )At] 

At^oo g[( Sl - S 2 )At] 

oo. (4.10) 



Since Eq. (14.101) contradicts the assumption that lim At^oo (s 2 At+t w ) lg\(s\— s 2 ) At] 

M =t m /p(At) 



is finite, the aging regime for experimental time scales does not exist in this case. 
Hence, if only this aging regime characterized by such g(t) exists, limAt->oo C(siAt + 

t w , s 2 At + t w ) is equal to or q. It implies that the work is independent of h(s). An 
example of such g(t) is g{t) = exp(t). 

Consequently, we conclude that if there are aging regimes characterized by g(r) of case 
2. or case 3. and the function C(/i) is not a constant and the FDT violation factor X is not 
equal to 0, the work for the slow process depends on the path of changing the external field 
h(s). In this case, contrary to the fundamental assumption of thermodynamics, the work for 
the slow process is not the difference of a state function. 

V. THE WORK IN THE AGING REGIME; RESULTS FOR 
SINGLE-CORRELATION-SCALE SYSTEMS AND 
MULTI-CORRELATION-SCALE SYSTEMS 

It has been explicitly checked on several disordered models that in the long-time limit 
two situations with different dynamical behavior seem to exist [l|]: 

• There are systems with only one correlation scale apart from the quasiequilibrium 
regime, which we call "single-correlation-scale systems". For these systems, the cor- 
relation function scales as C(t,t') = C[^(t)/^(t')] and the FDT violation factor X 
is a constant. Equilibrium states of these systems are solved by a one step replica 
symmetry breaking ansatz. Examples are the p-spin spherical spin-glasses || and a 
Lennard- Jones binary mixture which is a model of structural glasses in a glassy state 
||. The real spin-glasses such as AgMn seem to belong to this class since it is known 
that the correlation function is scaled with single scaling function £(t). 

• There are systems such as Sherrington-Kirkpatrick (S-K) model which have an infinite 
number of correlation-scales apart from the quasiequilibrium regime [[| , which we call 
"multi-correlation-scale systems". For these systems, ultrametricity in time holds for 
any correlation such that C(t,t') < q; C(t 1 ,t 3 ) = mm[C(t 1 ,t 2 ),C(t 2 ,t 3 )} when ti > 
t 2 > t 3 and £3 — ► 00. The FDT violation factor X is a nontrivial function of C. 
Equilibrium properties of these systems are solved by a full replica-symmetry breaking 
ansatz. 

In this section, we analyze the path-dependence of the work for the slow process in the 
aging regime for these two classes of systems by using the general results obtained in the 



previous section, Sec. [TV]. It is important to note that only these two classes of systems seem 
to exist PJ. 
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A. The single-correlation-scale systems 



For the single-correlation-scale systems, scaling of the correlation function C(t,t') 
C[£(t) I '£(£')] holds. Hence, the work for the slow process in the aging regime is given as 



W S /(AH) 2 = - f 1 d Sl h{ Sl ) r ds 2 h{s 2 ) 
Jo Jo 



xy{(7 hm -7 — -}}. (5.1) 

Three explicit choices of the scaling function £(t) have been proposed so far; £(£) = t,£(t) = 
exp[(£/r ) 1_a /(l — a)} and £(t) = exp[ln a (t/r )] |lj. We analyze below the properties of the 
work for the slow process for each three cases. 

1. When = t, the correlation function decays in power-law as found in the trap model 
pOU and the p-spin spherical spin-glasses . The slow-process limit of the correlation 
function in Eq. ( |5.1|) depends on si and s 2 as 

hm gfi^+O = S .l±£. (5.2) 



This is an example of Eq. (|4.5|) in previous section. 



2. When = expf^^;^) 1- "]^ < 1) [§TJ], the slow-process limit of the correlation 
function in Eq. ( |5.1| ) depends on Si and s 2 through the scaling function which behaves 
as 



r £{ Sl At + t w ) 



exp[^^(s! - s 2 )]. 



(5.3) 



This is an example of Eq. 



3. When = exp[ln a (t/ro)](a > 1) f25f , the slow-process limit of the correlation func- 



tion which appears in Eq. ( |5.1| ) depends on Si and s 2 through the scaling function 
which behaves as 

hm feAt + g =expM Sl - S2 )]. (5.4) 



This is another example of Eq. ( |4.8|) 

To my knowledge, only three types of the scaling function £(£) discussed above are known 
so far for single-correlation-scale systems. For each of three cases the slow-process limit of 
the correlation function in the integrand of Eq. (|4.1|) is a function of s% and s 2 . It means for 
the single-correlation-scale systems that the work for the slow process in the aging regime 
depends on the path of changing the external field h(s) and contrary to the assumption of 
thermodynamics the work for the slow process is not a difference of a state function. 
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In order to see more explicit form of the work for the slow process in the aging regime, 



we introduce the modified FDT Eq. (2J3) with constant violation factor X which has been 
verified to be valid for the single-correlation-scale systems M. 
From Eq. (|2.8|) , the susceptibility is given by 

x (t, t 1 ) = --±—C(t, £') + ^4tt7 ~ C ^~ 9 , (5-5) 
AV ' ' k B T e ff v ' ; k B T £ ff k B T ' v ' 



where C(0) = lim^oo C(t', t'). Hence, from Eq. (|4.1| ), the work for the slow process is given 
by 

H///Am2 C(0)-q 1 f 1 , dh{ Sl ) fi d/i(s 2 ) 

W B /(AH) = — — - - — - -tt / d Sl — -— I ds 2 



2k B T k B T e ff Jo ds\ Jo ds 2 
x{q- lim C[£( Sl At + tJ/£(s 2 At + t w )]}. (5.6) 



At— >oo 
t w ~At 



Since lim At->oo C(siAt + t w ,S2At + < g when si 7^ s 2 , the second term on the 

right hand side is non-zero. Hence, if T e -^ is finite the work for the slow process is not 
the difference of a state function since the work depends on the path of changing the field 
h(s). In addition, the value of the work is different from the value of the work for the slow 
process in the quasiequilibrium regime: the first term on the right hand side of Eq. ( p.6| ). 
In addition, it is also shown that since the correlation function is positive and smaller than 
the dynamical E-A order parameter in this regime, there are bounds for the value of the 
quasistatic work as 

C(0)-q q W a C(0)-g 

2k B T 2k B T e ff ~ (AH) 2 ~ 2k B T ' 1 ' ' 

Eq. ( |5.7| ) shows that the work W s in the aging regime is independent of the path of changing 
the field and coincides with the work W s in the quasiequilibrium regime when the effective 
temperature T e ^ is infinite, which holds for the 2-spin spherical spin-glass model |27j and 



a model of phase separation [11 



B. Possibility of observation of path-dependence 

In the previous subsection, we saw that the work for the slow process depends on the 
path of changing the external field when At ~ t w . We show here that the path-dependence 
of the work for the slow process is strong enough to be observed in experiments by evaluating 
the path-dependence of the 5-spin spherical spin-glass model. 

First, we describe the several pieces of information needed to evaluate Eq. (|5.6|) . It is 
known that the scaling of the correlation function, such that C(t,t') ~ q(t' /t) 1 and the 
modified FDT with the effective temperature, i.e. Eq. ( p.8|) , hold for the model ||. When 
k B T / J = 0.2 [PS|| , 7 ~ 0.25 has been obtained. The dynamical E-A parameter q and the 
effective temperature are given by q ~ 0.929 and T/T e *f ~ 0.227 when k B T / J = 0.2. In 
computation of these values, we uses the formulae: 
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FIG. 4. The relative difference between the works for the two paths (h(s) = s and h(s) = s 2 ) 
of changing the field infinitely slowly is plotted against the ratio of the waiting time to the mea- 
surement time for the 5-spin spherical spin-glass model. The path-dependence of the work for the 
slow process in the aging regime is strong enough to be observed in simulations for the model. In 
addition, we see that the path-dependence becomes larger, as one increases the measurement time 
with keeping the waiting time fixed. 
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In order to evaluate the magnitude of the path-dependence, we compute the relative 
difference between the works for two paths (h(s) = s and h(s) = s 2 ) which is defined as 
the difference between the works for the two paths divided by the average of the two works. 
The results obtained from Eq. (|5.6|) are as follows; When the ratio of the measurement time 
to the waiting time is 1 : 1, the relative difference is equal to 0.037, i.e. about 4%; When 
the ratio is 5 : 1, the relative difference is equal to 0.094, i.e. about 9%. 

Hence, we can say that the path-dependence of the work in the slow process can be 
observed for the 5-spin spherical spin-glass model. The relative difference is plotted against 
the ratio of the waiting time to the measurement time in Fig. |4]. It tells that contrary to 
intuition the path-dependence becomes larger, as one increases the measurement time with 
keeping the waiting time fixed. 



C. The multi-correlation-scale systems 

For the multi-correlation-scale systems such as S-K model, it is known that the ultra- 
metricity in time holds for all correlations such that C(t,t') < q in the long-time limit 

C(tt, t 3 ) = min[C(ti, t 2 ), C(t 2 , t 3 )\ when t l >t 2 > t 3 . (5.8) 
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By using this ultrametric relation Eq. ( |5.8| ), we prove below that the correlation function 
in the aging regime is a constant in the two cases (case 2 and case 3 in Sec. [TV]) where the 
work for the slow process can depend on the path of changing the external field h(s). Thus, 
for these two cases, the work for the slow process is independent of the path h(s). It has 
been proven in Sec. [TV] for the other two cases that the work for the slow process cannot 
depend on h(s). Thus, the work for the slow process for the multi-correlation-scale systems 
is independent of the path of changing the field in all the experimental time regimes. 



1. When g{r) = r, by substitution t\ — > t, t 2 — > fx't, t% — > fit in Eq. (|5.8| ) and assuming 
that fi' > fi, we see that 



Hence, 



lim C(t,fit) = minflim C(t,fj,'t), lim C(fi't, fit)]. (5.9) 

t — >oc t — *oo t — >oo 



lim C(fi't,fit) = lim Cifi't, —a't) = C(fiffi') (5.10) 



where C(fi) = lim t _ +0O C(t, fit). Thus, from Eq. (|5.iJ|), 

C(fi)=mm[C(fi'),C(fi/fi')]. (5.11) 

Since it implies that 

C"M = (S^' )if t f > " (5.12) 
I C (fjL ) otherwise, 

C(fi) is a constant and hence lim ^-*oo C(t + t', t') = C[u/ (1 + v)\ is a constant. 

2. When lim T _ >0 o'7~/5 , v r ) = and lim^oo g{r)/g{sT) = g r (s) (s < 1) exists, the ultramet- 
ric relation Eq. ( |5.8| ) is written as 



lim C(t + t' t') 

fi=t'/g( T ) 

min[ lim C(r + t',t'), lim C(ar + t',t')] (5.13) 



where a is a constant such that < a < 1. The second argument of min in Eq. ( |5.13| ) 
is rewritten as 

lim C(aT + t',t')= lim Ciar + t't') 

l*=t* fg(j) H=t f I g(ar I 'a) 

= lim C(T + t',t') 

fj,=t f fg{r/a) 

lim C(t + t',t') 

^=f/W«)sH] 



lim C(r + *',*') 

fj,g r (a)=t f /g(r) 

= C[fig r {a)\, (5.14) 
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where C(/x) = lim j" t 7*°° } C{r-\-t',t'). Similarly, the first argument of min in Eq. (|5.13|) 
is rewritten as 

lim C(t + if, ar + t') = lim C(r + t', ar + if) 

f*=t'/ 9 (T) M=(t'+ar)/s(r) 

= lim C[(l-a)r + t',t'] 
= C^l - a)] (5.15) 
where the last equality is shown from Eq. (|5.14j) . Consequently, 



C(//) = wia{C\pig r (l - a)],C\pig r (a)]}. (5.16) 
Since qJa) > 1, = lim t-+oo C(r + t' £') is a constant. 

A«=t'/s(r) 

We proved above only that the correlation function is a constant within a correlation- 
scale since ultrametricity in time is valid within a correlation-scale. On the other hand, 
even for the multi-correlation-scale systems, the value of the correlation function becomes 
when the time difference is much larger than the waiting time. So, value of the correlation 
function must depend on an infinite number of correlation-scales. Hence, there is possibility 
that there is an infinite number of experimental time regimes characterized by different 
values of the work for the slow process. 



VI. DISCUSSION AND CONCLUSIONS 

Summarizing our results, for systems with aging we have considered the experimental 
time scale dependence of the work exerted by modulating the external field in an isothermal 
process. We have shown that there are two experimental time regimes characterized by 
different behavior of the work. In the quasiequilibrium regime (At <C t w ), the work for 
the slow process is independent of the path of changing the external field and the value 
of the work for the slow process is different from the change of the Helmholtz free energy 
obtained by statistical mechanics with the assumption of ergodicity. When the process is 
not infinitely slow, the difference between the work for non-slow process and the work for 
the slow process is positive for any path of changing the filed, i.e. the principle of minimum 
work holds. In addition, we derived the dependence of the difference on the measurement 
time from the long-time behaviors of the correlation function. 



In the aging regime (At ~ t w ), whose precise meaning is given in Sec. |V|, contrary 
to a fundamental assumption of thermodynamics, the work exerted in the slow process on 
the single-correlation-scale systems for which the FDT violation factor X is not equal to 
depends on the path of changing the external field. It was examined that for the 5-spin 
spherical spin-glass model the magnitude of the path-dependence is strong enough to be 
observed in experiments. The conditions when the work for the slow process depends on the 
path is also derived generally. On the other hand, for the multi-correlation-scale systems, the 
work for the slow process is independent of the path of changing the field. However, the work 
for the slow process still may depend on history of a system since there is possibility that 
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there is an infinite number of experimental time regimes characterized by different values 
of the work for the slow process. In order to attack on this problem, information on time 
scale of the aging regime of multi-correlation-scale systems beyond time reparametrization 
invariance is necessary 
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It is important to note that these results have a much broader range of validity, because 
they were derived without any assumption on microscopic properties of a system. 

We consider the implications of our results on thermodynamics for systems with aging. 
In the quasiequilibrium regime, our results are the same as the consequences of thermo- 
dynamics, where the work for the slow process is independent of the path of changing the 
external field and the work for a non-slow process is larger than that for the slow process. 
Only discrepancy is difference between the value of the work for the slow process and that 
of the free energy obtained by statistical mechanics with the assumption of ergodicity. How- 
ever, this discrepancy does not mean that the framework of usual thermodynamics is invalid 
since the free energy can be defined by the work for the slow process. This suggests the 
usual framework of thermodynamics is valid in the quasiequilibrium regime even though the 
system is out of equilibrium. 

We considered only the properties of the work in the linear response regime. Thus, in 
order to prove validity of usual thermodynamics in the quasiequilibrium regime, the relation 
between the heat and the entropy has to be considered and extension of our analysis outside 
the linear response regime has to be performed. For these purposes, it would be interesting 
to analyze Langevin dynamics of spin-glass models with stochastic energetics |T9|j30| . 

Since the work for the slow process in the quasiequilibrium regime is independent of 
history of a system in linear response regime, it gives a fundamental state function for 
thermodynamics in the quasiequilibrium regime and is worth being measured in experiments. 
We discussed an experiment to measure it. Our results in Sec. fTJ for the measurement time 
dependence of the work shows how long measurement time is needed to reach the expected 
accuracy. 

However, since the waiting time can be long but should be finite in reality, the measure- 
ment time has to be finite in order to keep the experimental time scale in the quasiequilibrium 
regime. Otherwise the value of the work does not come close to that of the work for the 
slow process in the quasiequilibrium regime since the experimental time scale is in the aging 
regime. Hence, one cannot expect that the value of the measured work coincides with that 
of the work for the slow process with arbitrary accuracy. In other words, the accuracy of 
measurement of the state function is bounded by the length of the waiting time. 

On the other hand, in the aging regime, our results imply that in contrast to the usual 
framework of thermodynamics it is essential to consider the experimental time scales such 
as the waiting time and the measurement time. Furthermore, history of a system as a state 
variable is necessary to describe thermodynamic properties of glassy materials. Although 
constructing thermodynamics for the aging regime is a tough work, it is a fruitful and 
challenging problem to be investigated since some universal properties in the aging regime 
are known ||. Recently, a framework of thermodynamics for the aging regime when the 
temperature changes in time is proposed ||31| . The framework avoids treating the problem 
of path-dependence by fixing a path to change the temperature. It is not clear whether 
thermodynamics for a single path to change the temperature is useful. Its validity and 
usefulness should be checked. 
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APPENDIX A: DERIVATION OF THE LONG MEASUREMENT TIME 

BEHAVIOR OF W-W s 

In this appendix, we derive the long measurement time behavior of the deviation from 
the work for the slow process for four types of relaxation of the correlation function described 
in Sec. [DTI]. 

1. The case where lim^oo t[C{t) — q] exists 

In this section, we treat the case where lim^oo t[C(t) — q] exists, which includes the two 
cases where lim^oo t[C(t) — q] is equal to and where the limit is equal to a finite value c. 
We rewrite the deviation given by Eq. ( ft. ICQ by introducing K(t) defined as 



K(t)= f dr[C{r) - q]. (Al) 
Jo 

The deviation is reduced to 

+ J* 1 ds 2 ti'( Sl -s 2 )K(s 2 At)]. (A2) 

At first, we consider the deviation when lim^oo t[C(t) — q] is equal to 0. Here, in order 
to prove that the deviation is proportional to the inverse of the measurement time At, we 
show that the proportional constant is finite. The proportional constant is given as 

lim At-^—(W-W qs )= lim /* d Sl ^^[ti (ti)K(siAt) 

+ f * ds 2 h"( Sl - s 2 )K(s 2 At)}. (A3) 
Jo 

Since lim^oo K(t) exists in the case, K (t) is finite. Hence, the order of the limit At — > oo 
and the integration can be changed. Thus, the proportional constant is given by 



/ ds 


dh(s) 


Jo 


ds 



(A4) 



where is defined by Eq. ( |3.13|) . Consequently, it is proved that the proportional constant 
is finite and the deviation is proportional to 1/ At. 
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Next, we consider the deviation when lim^oo t[C(t) — q] is equal to a finite value c. In 
order to prove that the deviation is proportional to log At/ At, we show that the proportional 
constant is finite. Dividing the region of integration of Eq. ( |A2| ) into five regions, the 
proportional constant is given by 



lim 



At k B T{W -W t 



qs ) 



At-^oo log At (AH) 2 
lim 7-^r-MO) P ds l h'(s 1 )K(s l At) 

At^oo log At JO 



+ ti(0) I ds 1 h'{s 1 )K(s 1 At) 

+ f^ds 1 h'(s 1 ) r ds 2 h"{s x -s 2 )K{s 2 At) 
Jo Jo 

+ / d Sl h'( Sl ) [ At ds 2 h"(s 1 -s 2 )K(s 2 At) 
J ir t Jo 

+ I d Sl h'{ Sl ) ds 2 h"{ Sl - s 2 )K{s 2 At)\. (A5) 

"'St •''St 

We discuss the behavior of the function K(sAt) in each region. When s < 1 /At, from 
the definition of K{t), Eq. (|3.13| ), it is shown that 

KisAt) < f 1 dr\C(r) - q]. (A6) 
Jo 

Thus, K(sAt) is finite even when At — > oo. On the other hand, when s > 1/At, we express 
the function K(sAt) by using a function defined as C(t) — q = c/t + f{t),t > 1: 



KisAt) = [ SAt drlCM - q] 
Jo 

= [ 1 dr[C(r)-q]+ [ SAt dr[- + f(r)} 

JO Jl T 

= J dr[C(r) -q] + clogs + clog At + J' drf(r). (A7) 

Since lim^oo tf(t) = 0, the fourth term on the right hand side is finite even when At is 
infinite. Thus, the behavior of the function is given as K(sAt) ~ clog At when s > 1/At. 
By evaluating the each term of Eq. ( |A5| ) with using the above behavior of the function 
-fT(sAt), it is straightforward to show that 

At k B T /ttt xx r x f 1 j >dh(s). 2 . . . 

hm — . A TT ,J W - W ga ) = c / ds — — . (A8) 

At^oologAt (AHf y qJ Jo 1 ds 1 y J 

Consequently, it is proved that the proportional constant is finite. 
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2. The case where lim^tx, t[C(t) — q] does not exist 

In this section, we treat the case where lim^oo t[C(t) — q] does not exist, which includes 
the two types of long-time behavior of the correlation function, such that C (r) ~ q+c/ r a , < 
a < 1 and C(r) ~ q + cj logr. 

In order to find out the leading term in the behavior for the long measurement time, we 
divide the region of integration of Eq. ( 3.10Q into three regions as 

/■Si AAt)~ a rsi rl rsi 

ds\ I ds 2 = / dsi / ds 2 + / ds\ \ ds 2 



l(At)- a Js 1 -(At)- a 
rl r s!-(At)- a 

+ / dsx / ds 2 , (A9) 

J(At)- a JO 

where we assume that the parameter a satisfies < a < 1. Since all the integrands (the 
time derivative of the field and the correlation function) are finite, it can be shown that the 
first two terms on the right hand side are bounded as 

| the first term| < const x (At)~ 2a 

and 

| the second term| < const x (At)~ a [l - (At) - *]. 

In addition, we notice that since s 2 < si — (At)~ a in the region of integration of the third 
term, (si — s 2 )At is large when the measurement time At is large. Thus, the third term is 
determined by the long-time behavior of the correlation function. 

At first, we treat the case where the correlation function behaves at a long time as 
C(t) ~ q + c/r~ a , < a < 1. We divide the region of integration of the third term of Eq. 
|) into three regions as 

dsi / ds 2 = / ds\ I ds 2 — / ds\ I ds 2 

(At)~ a Jo Jo Jo Jo Jo 

dsi r ds 2 . (A10) 

(A)- a J sl -(At)- a 

Since all the integrands are finite, it is shown that the last two terms are bounded by 

| the second term| < const x (At)~ a(2 ~" )-a 

and 

| the third term| < const x (At)" a " o(1 " a) . 

On the other hand, by using the long-time behavior of the correlation function, the first 
term is given as 

/ygfe) r^ dh } S2) , c w A V < ai d 

Jo da x Jo ds 2 (si - s 2 ) a (At) a y ' 
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Thus, we find that this term is 0[(At)~ a }. Assuming the free parameter a is larger than the 
exponent a, one sees that the leading term is Eq. ( |A11| ) by comparing the orders of all the 
terms in Eqs. ( |A9| ) and ( |A10| ). Hence, the long measurement time behavior of the deviation 



is given by Eq. ( |3.15|) . 



Next, we consider the case where the long-time behavior of the correlation function is 
given by C(r) ~ q + c/ logr. Again, we start by evaluating the third term on the right hand 
side of Eq. ( |A9| ) . By using the long-time behavior of the correlation function, the term is 
given by 

/' dsi d Jpl r-^- ds /Jp>_ " (A12) 

J(At)- a dsi Jo ds 2 log(si — s 2 ) + log At 

Since in the region of integration the following inequality holds: 

logOi - s 2 ) 



log At 



< a, 



by neglecting numbers of Oiaj log At) the third term of the right hand side of Eq. (|A9|) is 
reduced to 

/' ds /Jpl r~^ ds /jpic (A13) 

J(At)- a dsi Jo ds 2 log At 

Furthermore, by neglecting numbers of 0[(At) _a / log At], the term is reduced to 

1 , dh(si) , dh(so) c c 

dsi , ' / ds 2 ; 2J — = — — . (A14) 

o dsi Jo ds 2 log At 2 log At 

Consequently, one sees that the leading term is the above one by comparing the orders 
of all the terms in Eqs. (|A9|) and ( |A10|) . Thus, the long measurement time behavior of the 
deviation is given by Eq. fl3.16| ). 



Finally, in order to determine the order of the neglected term we determine the value 
of the free parameter a, since the neglected term contain the free parameter a. Since in 
terms of a the order of the neglected terms are given as Oiaj log At) and 0[(At)~ a ], we 
can determine the value of a so that the value of the neglected terms are minimized. Thus, 
assuming that the neglected term is given by aa/logAt + /3(At)~ a where a and (3 are 
constants, the value of a is determined by 

^ a i^ + " (A ^i = °- (A15) 

Since the solution is a ~ 21og(log At)/ log At, one sees that the order of neglected term is 
0[log(log At)/(log At) 2 ] as shown in Eq. (13T6Q . 
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